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Numerical simulations of variants of the CDM model with different cosmological 
parameters are used to compare statistical measures such as mass spectra, merger 
processes, and autocorrelation functions, for different models with relevant observa- 
tions. The degree of virialisation of the halos is checked, and also which properties 
distinguish recent mergers. Mergers occur mostly in deep potential wells and mark 
the most important structure formation processes. As consequence, the autocorre- 
lation function of merged halos has a higher amplitude and is steeper than that of 
the virialized clusters. This effect can also be seen for ultraluminous IRAS galaxies 
which are thought to be results from ongoing merging events. 

1 Introduction 

Gravitational instability is accepted to be the basic driving agent of structure 
formation on large scales. Combined with the CDM model it leads to the pic- 
ture of hierarchical clustering with wide support from deep galaxy and cluster 
observations. The most massive virialized objects we observe in the universe 
are clusters of galaxies. We can derive lots of cosmological information by 
investigting these objects in detail: Since they are just in the process of rising 
from rare high peaks in the primordial density field, non-linear effects have 
not erased all the information about their formation. Additionally the mass 
and abundance of galaxy clusters gives us a hint what the amplitude of the 
fluctuations in the primordial density field had been. 

In order to derive reliable physical properties of galaxy clusters from nu- 
merical simulations one has to use large simulation boxes and a high resolution 
TV-body code, because the overall large scale structure is connected to the sub- 
structure of halos (Colberg et al.l997). 

2 Cosmological models 

We base our analysis on three cosmological models all with an amplitude as 
determined by the 4- year COBE experiment, cp. Bunn & White (1997). First 
we take the standard CDM-model, with critical mass density JIq = 1 and a 
dimensionless Hubble constant h = 0.5, which is used as a reference in spite 
of the difficulty in reproducing both large scale and small scale clustering. 
More realistic variants of the CDM model have a lower matter content, as 
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Table 1: Physical properties of the numerical simulations. The box size L is given in /i ^Mpc, 
the particle mass nip in units of 10^^ /i~^M0. 
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ACDM which has flo = 0.3, a cosmological constant of Xq = 0.7, and h = 0.7. 
For this model the shape parameter of the power spectrum F = rtgh — 0.21 
better fits the constraints from galaxy and cluster clustering, (e.g. Peacock 
& Dodds 1996, Einasto et al.l998). The models were normalized to the full 
4-year COBE anisotropy using the Boltzmann code CMBFAST developed by 
Seljak & Zaldarriaga (1996) and assuming a baryon content of = 0.0125 
suggested by big bang nucleosynthesis. The mass variance of the linear input 
spectrum at a scale of 8/i^^Mpc is given by erg- Finally we used a less extreme 
open model OCDM with JIq ~ 0.5 which promises realistic large scale galaxy 
clustering, and which has the same mass variance at 8/i~^Mpc as galaxies. All 
information about the models can be found in Table ^ and the corresponding 
power spectra are plotted in Fig |l|. 
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Figure 1: Power spectra for the models specified in Table 
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3 Numerical Techniques 



3.1 N-body simulations 

To set up the numerical simulation of a gravitating system of dark matter 
particles one has to impose density fluctuations onto a distribution that repre- 
sents a homogeneous and isotropic universe. Therefore a glass-like distribution 
is used for producing a starting realization (White 1993). To give an impres- 
sion of the improvement using glass starting places as input for the following 
numerical simulation I have plotted in Fig |l| a slide through the initial particle 
distributions, using grid starting places and glass starting places, respectively. 
It can easily be seen that the input power spectrum is much better represented 
using the glass data. 
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Figure 2: Representation of the inital power spectrum (SCDM model) using grid starting 
places (left panel) and glass starting places (right panel), respectively. 

The evolution of the initial particle distribution is simulated using a modi- 
fied version of the adaptive P'^M code (Couchman 1991) which is able to allow 
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Table 2: Linking lengths II for the different models used with the friends-of-friends algorithm 
and the percentage of identified unvirialized particle groups. 
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for background expansions other than the Einstein-de-Sitter universe. The 
adaptive part of the code counter balances the dramatic slowdown caused by 
direct summation over nearby particles by placing refinement grids at the high 
density regions as clustering evolves. The force calculation for particles be- 
longing to these refinements is performed by the same P'^M code on a finer 
grid with isolated boundary conditions. 

3.2 Identification of galaxy clusters 

Groups of galaxies and galaxy clusters are identified with groups of dark 
matter particles defined with a standard friends-of-friends algorithm (Davis 
et al.l985). This methods collects particles in groups whose spatial separa- 
tion ru is smaller than II (Zinking Zength) times the mean inter particle spac- 
ing Aa;: 

Til < II ■ Ax. (1) 

At this point the question of the appropriate value for the linking length II 
arises. This parameter should be chosen in order to identify gravitational 
bound objects in the simulations. If we assume spherical halos with an isother- 
mal density profile p(r) cx 1/r^ and the usual value for the linking length U — 
0.2 we will find objects with an overdensity of S = {p — 'p) /'p ^ 180 or higher 
(Lacey & Cole 1994). This value corresponds to 5th ~ 178 calculated for 
virialized objects under the assumption of a spherically symmetric collapse in 
an Einstein-de Sitter universe. In order to get the correct linking lengths for 
gravitationally bound objects in open and A-universes, we use the following 
formula: 

/ 170 \ 1/3 

'^(^™) = 0.2 ^ ^ . (2) 

The values of II listed in Table || are based on eq. || where the values of 
(5th(^^07 ^^a,o) are calculated from the formula found in Katayama (1996). 
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However, even though the theory tells us that we should get virialized 
objects with a mean overdensity as given above, we cannot be sure that the 
particle groups we find by using the linking lengths from Table |^ in our friends- 
of- friends scheme are really physically virialized objects. In order to check this 
I have explicitly tested the virial theorem |-Epot| = 2 -Ekin for each individual 
halo. 




Figure 3: Virial theorem for clusters in the ACDM model (units are arbitrary). 

Fig. ^ shows the virial relation as a scatter plot for all particle groups 
identified in the ACDM model at a redshift of z = 0. To emphasize the impor- 
tance of high (spatial) resolution we have calculated the same model with a 
PM code that was run with exactly the same starting configuration as used for 
the AP^M code. The velocity dispersions in the PM simulation are throughout 
underestimated and lead to a correlation of the order |i?pot| ~ 8.2 i?kin- The 
'lower' virial relation found in the AP"^M simulation, |i?pot| ~ 1-8 i?kin, can be 
understood as an influence of an outer pressure of radially infalling particles 
into the halos (Cole & Lacey 1996). The obvious bend of the virial relation at 
low energies towards smaller -Bkin proves that the surface influence decreases for 
groups of small particle numbers and the expected behaviour l-Epotl ~ 2 -Ekin 
is reached. An opposite effect concerns a much smaller number of groups with 
velocity dispersions which are too high compared to their potential energy. 
These clusters are considered unvirialized and gave rise to question what ef- 
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fects are responsible for this deviation from virial equilibrium. Therefore they 
are treated separately in the following analysis, and are investigated in de- 
tail. Table ^ also shows the percentage of unvirialized objects among all FOF 
clusters for the different cosmological models. It becomes obvious that the 
percentage of unvirialized groups is highest in SCDM and lowest in OCDM. 
This is an representation of the different absolute abundance evolution in the 
different models. 



4 Properties of (un)virialized halos 

One of the most basic properties of a galaxy cluster is its mass M . Therefore 
I would like to present the cumulative massfunction dn/dM in comparision to 
the prediction from the Press-Schechter theory (Press & Schechter 1974) and 
to observational data compiled by Bahcall & Cen (1993). 




Fig. ^ shows that the ACDM model best fits the observation and is very 
close to the prediction made by Press & Schechter. The PS curves for the 
SCDM and OCDM are left out for clarity; they match the numerical as well 
as the ACDM model. 

In the following the nature of the unvirialized halos will be investigated in 
more detail. 
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First of all I will discuss the consequences of merger events for the shape 
of the resulting groups. To this aim the eigenvalues of the inertia tensor for 
each FOF group are calculated, which are given in the order a > b > c. In 
Fig. H we show a scatter plot of the ratios of the eigenvalues for the ACDM- 
simulation at z = 0. In this plot, spherical groups are situated near the origin, 
oblate clusters in the upper part and prolate clusters in the lower right corner. 
It is well know (Dubinsky 1992, Warren et al.l992, Lacey & Cole 1994) that 
hierarchical clustering leads to triaxial ellipsoids with a typical axis ratio of 
6:4:3. The unvirialized halos inhabit more the right part of the diagram, again 
characterizing the soft merging, i.e. the elongation of the groups due to tidal 
interaction of the progenitors which marks the direction of the encounters, 
and in some cases to an elongation of the clusters in a second direction due 
to non-central encounters. Typically unvirialized clusters have an axis ratio of 
8:4:3. 
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Figure 5: Axis ratio of virialized (dots) and unvirialized groups (diamonds) for the ACDM- 

simulation at 2 = 0. 



Next I concentrate on the correlation of velocity dispersion and the mass 
of virialized and unbound groups. A tight correlation of both quantities for 
virialized objects is expected if we assume an isothermal sphere for the halo 
and a cutoff at constant density, which can easily be taken to be the mean 
cosmic density of the halos at the time when they form or suffer their last big 
merger. Since most halos form quite recently this cutoff is almost constant. 
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and we get 

(Tv = K/\/2 cx Mi (3) 

As can be seen from Fig. ^ this is represented very well in our simulations 
with a large scatter for the light clusters. Fig. ^ also suggests that the unbound 
halos mostly lie at the low mass end of the distribution with high inner velocity 
dispersion. Cole & Lacey (1996) argue that the tail of groups with cr„ > VI3/V2 
are likely objects that belong to larger virialized structures. This may explain 
why we identify them as unvirialized objects. But if we repeat our analysis for 
objects obtained using a greater linking length II we derive the same results. 
Even for FOF groups identified with a linking length corresponding to an 
overdensity S ^ (5th we find unbound groups in the tail ay > Vc/V2. 
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Figure 6: Relation between velocity dispersion and the mass for the ACDM model. The 
virialized halos are marked with points and the objects that are not in virial equilibrium are 
marked using diamonds. The solid line is a fit to the scaling relation cTv Ms . 



More interesting than these low mass groups are the small number of halos 
lying beneath the virial relation = Vc/V^. They represent groups in the 
process of merging which are unbound because they are too extended. In most 
cases they represent halos with more then one centre which are connected by 
slight tidal bridges and similar structures. These systems are not very frequent, 
i.e. they do not survive for a long time, but they mark most interesting places 
in the simulation box. And it can be expected that they are sites for active 
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structure formation processes in nature. Fig. g (at the end of this article) shows 
the most massive unviriahzed cluster in the ACDM simulation at z = 0.1 and 
its two almost equal mass progenitors at z = 0.2. Both these progenitors are 
virialized objects just like the final merger product at a redshift of z = 0. 

5 Correlation functions 

One of the basic constraints of cosmological models is the shape and the am- 
plitude of the two-point correlation function. For many years it has been the 
standard way to describe the clustering of galaxies and galaxy clusters. The 
assumption that galaxies (and clusters) only form from high-density regions 
above some threshold value Sc leads to a correlation of points exceeding this 
value 6c that is enhanced in comparison to the dark matter correlation function 
(Kaiser 1984). 




Figure 7: Correlation functions for the ACDM model. Histograms from below (at r > 
3/i~^Mpc) denote dark matter particles, virialized clusters and unvirialized groups, both 
with a mass cut corresponding to 50 particles. The dashed lines are fits with parameters 

specified in Table ^ 

In Fig. 1^ we show the correlation functions for the ACDM model. We 
find that the correlation between virialized groups containing more than 50 
particles corresponding to a mass cut of 4 • 10"'^'^ h~^MQ lies by a factor two 
over the correlation function of dark matter. We fitted standard power laws 
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Table 3: Fit parameter of the correlation functions for dark matter (first column), bound 
clusters (second column) and unbound clusters (third column) for three cosmological models. 
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e(r) = (ro/rr (4) 

which are shown in the Fig. |^ as dashed hues. Galaxy groups in the ACDM 
model have a correlation length ~ 9 h~^Mpc and a slope 7 — 1.7. The cor- 
relation function turns negative beyond 30 h~^Mpc. Fits for the other models 
are given m Table |. To get comparable values for the different cosmologies 
we fixed the cluster number density to n — 10"'^/i'^Mpc~'^ and used only the 
N = nV most massive clusters for the calculation of the correlation function. 
Obviously the amphtude in SCDM and OCDM are smaller than in the ACDM 
model, i.e., these models have too small power on large scales for describing 
the clustering of groups and clusters of galaxies. 

The correlation function of the unvirialized groups is also of interest. They 
are much more strongly clustered than the virialized objects. Also the slope of 
the correlation function is very steep, 7 = 2.5, independent of the cosmological 
model, cp. Table ^. This verifies that mergers are occurring at particular places 
in the universe, and that these processes are highly correlated. A similar result 
has been found for ultraluminous IRAS galaxies (Gao 1993). 

6 Conclusions 

We have collected a number of reasons for looking at virialization as a rea- 
sonable criterion for identifying groups and clusters of galaxies in numerical 
simulations. The virial theorem well described an overwhelming majority of 
the halos in our simulations. A slight bend in the virial relation can be as- 
cribed to the effective pressure of the permanent spherical accretion stream in 
the halos. Special attention has been devoted to the unvirialized halos. We 
have shown that they are characterized by quite recent soft mergers which 
lead to more anisotropic halos and are strongly correlated over scales of up 
to 40 /i^^Mpc. In our simulations we find a very fast virialization of the ha- 
los which leads to the small percentage of such objects at any given time. It 
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must be checked whether the viriaHsation time depends on the relatively low 
mass resolution in our big simulation boxes. The stronger anisotropy of the 
merger products leads to typical triaxial ellipsoids with axis ratios of 8:4:3. 
The merging processes are characterized by basically central encounters. Not 
much angular momentum is transferred to the merger product which leads to 
a self-similar growth of the rotation of the halos, (e.g. Knebe 1998). 

The large differences in the correlation functions of the models makes it 
worthwhile to further compare differences in the large-scale matter distribu- 
tions of the different models. 
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Figure 8: Projections of a cube of (lO/i^'^Mpc)^ of the ACDM simulation box at redshifts 
z = 0.2, z = 0.1 and 2 = 0, containing the most massive unvirialized cluster. Only the 
particles in the cluster and in the two most massive progenitors are shown which are both 

virialized objects. 
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